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Abstract. We evaluate the matrix elements {r^) for the n-dimensional harmonic oscillator in 
terms of the dual Hahn polynomials and derive a corresponding three-term recurrence relation and 
a Pasternack-type reflection relation. A short review of similar results for nonrelativistic hydrogen 
atom is also given. 



1. An Introduction 



The purpose of this note is to present a simple evaluation of the expectation values (r^) for the n- 
dimensional harmonic oscillator in terms of the dual Hahn polynomials by direct integration. Other 
methods of solving similar problems in elementary quantum mechanics appeal to general principles 
and involve the Hellmann-Feynman theorem (see |6j, [7], [S] and references given there, textbooks 
EZ], [22], [m, HO]), commutation relation ([23], [SD], and dynamical groups ([3], |1], [S], 

[TU] . [TB] . [37], [S]). Our approach allows to study these expectation values with the help of the 
advanced theory of classical polynomials [2] , [25] , and [35] . We recall also that the first-order of the 
time-independent perturbation theory equates the energy correction to the expectation value of the 
perturbing potential. Thus expressions of the form (r^) gain utmost importance. 

The plan of the paper is as follows. We review results on expectation values (r^) for the nonrela- 
tivistic hydrogenlike wave functions first and then extend them, in a similar fashion, to the case of 
n-dimensional harmonic oscillator. An attempt to collect the available literature is made. 



The problem of evaluation of matrix elements (r^) between nonrelativistic bound-state hydro- 
genlike wave functions has a long history in quantum mechanics. An incomplete list of references 

includes [1], m, n, m, m, m, m, m, m, m, m, m, m, m, m, m, m, m, m, 

[10], [IS], [H], [50], [51], and [56] and references therein. Although different methods were used in 
order to evaluate these matrix elements, one of possible forms of the answer seems has been missing 
until recently. In Ref. [39] the mean values for states of definite energy 
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2. Expectation Values (r^) for Coulomb Problems 
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in the nonrelativistic Coulomb problem have been evaluated in terms of the Chebyshev polynomials 
of a discrete variable tk {x, N) = hf' {x, N) originally introduced in Refs. [51] and [52] . Their 
extensions, the so-called Hahn polynomials, introduced also by P. L. Chebyshev [53] and given by 

"t'H-,N) = (-D'LW^ (2.2) 



X 3F2 



—k, a + P + k + 1, —X 
1-N 



were rediscovered and generalized in the late 1940s by W. Hahn. (We use the standard definition 
of the generalized hypergeometric series throughout the paper [3], [2T].) 

The end results have the following closed forms 

(r'^-^> = ^(^)'"4(--/-l,-2/-l), (2.3) 



when k = 0,1, 2, ... and 



2Z ^ ^"^^ 



k+2/ n \ I tkin-l -1,-21-1), (2.4) 
j.K+z ^ 2n \nao J 

when k = 0,1,..., 21. Here = fi^/me^ is the Bohr radius. Equations fl2.ip - fl2.3l) reflect the 
positivity of the matrix elements under consideration [19]. 



The ease of handling of these matrix elements for the discrete levels is greatly increased if use is 
made of the known properties of these classical polynomials of Chebyshev [21], [25], [35], [36] and 
[51] . [52], [53]. The direct consequences of these relations are an inversion relation: 



) (2.5) 



,r^+V \nao) {2l + k + l)\ ^ 
with < A; < 2/ and the three-term recurrence relation: 

<'■'> - ^^1^ (If) - (H)^ (^'-^> (^-f') 

with the initial conditions 

-) = ^, 1 =1, 2.7 
which is convenient for evaluation of the mean values (r^) for k >1. 

The recurrence relation (12. 6p was originally found by Kramers and Pasternack in the late 1930s 
[27], [38], and [39]. The inversion relation (12.51) . which is also due to Pasternack, has been redis- 
covered many years later [15], [33] (see also [22] and [23] for historical comments). Generalizations 
of (I2.5p - (I2.6I) for off-diagonal matrix elements are discussed in Refs. [H], [37], [23], [15], and [50]. 
The properties of the hydrogenlike radial matrix elements are considered from a group-theoretical 
viewpoint in Refs. [3], [IS], and [31]. Extensions to the relativistic case are given in [T], [T7] . 
and [IHI (see also references therein). 



In a retrospect, Pasternack's papers [M], [SH] had paved the road to the discovery of the continuous 
Hahn polynomials in the mid 1980s (see [26], [19] and references therein). 
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3. Evaluation of (r^) for Harmonic Oscillator in EJ^ 
The stationary Schrodinger equation for the n-dimensional harmonic oscillator 

m=E% i/ (3.1) 

s = l ^ A / 

can be solved by separation of the variables in hyperspherical coordinates. The normalized wave 
functions have the form 

m [x) = ^>MKu (r, 9) = Yku {^) Rnk (r) , (3.2) 
where Yxu {^) are the hyperspherical harmonics associated with a binary tree T, the integer number 
K corresponds to the constant of separation of the variables at the root node of T and v = 
{Zi, Z2, ••• , Zp} is the set of all other subscripts corresponding to the remaining vertexes of the binary 
tree T (see [35], [IB], [1^, and fS3] for a graphical approach of Vilenkin, Kuznetsov and Smorodinskii 
to the theory of spherical harmonics). 

The radial functions are given by 



where (^) are the Laguerre polynomials [36] and the corresponding energy levels are equal to 

E = EN = N + n/2, {N -K) /2 = k = 0,1,2,... . (3.4) 

(See also Refs. [30], [31], [35], [H], and [16] for group theoretical properties of the n-dimensional 
harmonic oscillator wave functions.) 

The expectation values under consideration 

= [ ^%KAr,Q)rP^^Ku{r,n) dv (3.5) 

POO 

Jo 

we use orthonormality of the hyperspherical harmonics, are derived with the aid of the integral |37j . 
TSI: 

POO 

Jls = / e-^x"+^ (x) Li (x) dx (3.6) 
Jo 

T{a + s + l)T{l3 + m + l)V{s + l) 



m\ {n — m)\ F (/5 + 1) F (s — n + m + 1) 



— m, s + P — a — s 
X 3^2 I ^ , , I , n > m, 

p + 1, n — m + 1 

where parameter s may take some integer values. Similar integrals have been discussed in Refs. [1] , 
[TT] , [H] , and [2B] ■ A simple evaluation of this integral is given in Ref . [IS] . 

Indeed, 

, p. _ {(N-K) /2)! l2+K+nl2-l ( jK+n/2-1 

^'^^-F((iV + ii' + n)/2) X ^ (^-7) 
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where we replace r = ^. As a result, 

r(/f + (n + p)/2) / (A- - A') /2. P/2 + L -P/2 . \ 
^'^ ' " r (K + n/2) A- + n/2, 1 ' ^ j ■ (^■'*' 

provided that p + n + 2K > 0. Connection with the dual Hahn polynomials given by ^35] 

m! 

— m, a — s, a + s + 1 

X 3^2 I 1 + ^ _ 5 1 + a + c ' ^ 



r(i\:+ (n + p)/2) (0) fP fP , A n ^ ^ 



is as follows 

Thus the direct calculation of (r^) for each admissible p involves a well-known special function 
originally introduced in Ref. p^. Basic properties of the dual Hahn polynomials are discussed in 
[35] . We shall use one of them in the next section. See also Refs. |T6], [4lj, and [4j for group- 
theoretical methods of evaluation of similar matrix elements. We shall elaborate on the group- 
theoretical meaning of Eq. (13.81) later. 

4. Three Term Recurrence Relation 
The difference equation for the dual Hahn polynomials has the form 

where A/ (s) = V/ (s + 1) = / (s + 1) - / (s) , a; (s) = s (s + 1) , xi (s) = x {s + 1/2) , and 

a{s) = {s~a){s + b){s-c), (4.2) 
a(s) + r(s) Vxi(s) =(t(-s-1) (4.3) 

= (a + s + 1) (6 - s - 1) (c + s + 1) , 
Xm = m. (4.4) 
It can be rewritten as a recurrence relation 

a (-S - 1) Vx (s) y{s + l) + a (s) Aa; (s) y (s - 1) (4.5) 
+ (A„Ax (s) Vx (s) Vxi (s) - cr (-S - 1) Vx (s) - a (s) Ax (s)) ?/ (s) = 0. 
See [24], [25] and [35] for more details on the properties of the dual Hahn polynomials. 



Equations fl3.10p and (14.51) results in the three-term recurrence relation 

(p + 2)(r^'+') = {p + I) {2N + n) {rP) (4.6) 

+P (^^^ _ (ir _ 1) (ir + p - 1)^ 
for the matrix elements. Such recurrence relation requires two initial results 
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in order to generate all expectation values. 
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The cases n = 1, 2, 3 are well-known and usually discussed with the help of the hypervirial theo- 
rems (see, for example, Pj, [H], [20], [7], [ID], and references therein.) Our approach emphasizes the 
simple fact that the matrix elements recurrence relations (12. 6p and (14.61) for the hydrogen atom and 
for the n-dimensional harmonic oscillator, respectively, have the same "special functions" nature, 
namely, they occur due to the difference equation for the corresponding dual Hahn polynomials. 
We leave further details to the reader. 



5. A Pastern ACK- Type Inversion Property 
From equation (13. 8p one gets 

(r--) - r(^+(^-p)/2-l) 

^~ T{K+{n + p)/2) ^^■^> 

(for all the convergent integrals) in the spirit of Pasternack's inversion relation (12.51) valid in the case 
of the hydrogen atom. It simply reflects an obvious fact that the argument x (p) = {p/2 + 1) p/2 
of the dual Hahn polynomials in ( 13.101) obeys a natural symmetry x {—p — 2) = x (p) on the corre- 
sponding quadratic grid. We were unable to find relations ( 13.81) . ( 13.101) . and ( 15. ip in the available 
literature. 



6. A Conclusion 

Special functions appear in numerous applications of mathematical and physical sciences and 
their basic knowledge is a necessity for any theoretical physicist. As Dick Askey mentioned once in 
his talk on mathematical education, 'We should not lie to our students, but we do not have to tell 
them the whole truth." The same is true about teaching elementary topics in quantum mechanics. 
As we can conclude from this short note, even if "one is able to evaluate average values by appeal to 
general principles rather that by direct integration [20]", important relations with the well-known 
special functions will be out of the picture due to the limiting mathematical tools. Every instructor 
has to find a way how to resolve this contradiction. 
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